%QEEG YBOCS for responsive & None responsive
% Delta
% Post
clc
clear all
load chanlocs19

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  Pre Reponsive
% load power_Deltaband_d219Preeo 
% Pxx_PreNone1=Pxx6;
load power_Deltaband_d312Preeo 
Pxx_PreNone2=Pxx6;
load power_Deltaband_d313Preeo 
Pxx_PreNone3=Pxx6;
load power_Deltaband_d112Preeo 
Pxx_PreNone4=Pxx6;
% load power_Deltaband_d316Preeo 
% Pxx_PreNone5=Pxx6;
load power_Deltaband_d326Preeo 
Pxx_PreNone6=Pxx6;
load power_Deltaband_d115Preeo 
Pxx_PreNone7=Pxx6;
% load power_Deltaband_d116Preeo 
% Pxx_PreNone8=Pxx6;
load power_Deltaband_d117Preeo 
Pxx_PreNone9=Pxx6;
% load power_Deltaband_d118Preeo 
% Pxx_PreNone10=Pxx6;
load power_Deltaband_d110Preeo 
Pxx_PreNone11=Pxx6;
% load power_Deltaband_d120Preeo 
% Pxx_PreNone12=Pxx6;

Pxx_PreNone_Delta1=[Pxx_PreNone2;Pxx_PreNone3;Pxx_PreNone4;Pxx_PreNone6;Pxx_PreNone7;Pxx_PreNone9;Pxx_PreNone11];
Pxx36=mean(Pxx_PreNone_Delta1);

subplot(1,3,1)
topoplot(Pxx36,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Responsive-Delta Band';'Before'})
h=colorbar;
h.Location = 'southoutside';

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  Pre Nonreponsive
load power_Deltaband_d219Preeo 
Pxx_PreNone1=Pxx6;
% load power_Deltaband_d312Preeo 
% Pxx_PreNone2=Pxx6;
% load power_Deltaband_d313Preeo 
% Pxx_PreNone3=Pxx6;
% load power_Deltaband_d112Preeo 
% Pxx_PreNone4=Pxx6;
load power_Deltaband_d316Preeo 
Pxx_PreNone5=Pxx6;
% load power_Deltaband_d326Preeo 
% Pxx_PreNone6=Pxx6;
% load power_Deltaband_d115Preeo 
% Pxx_PreNone7=Pxx6;
load power_Deltaband_d116Preeo 
Pxx_PreNone8=Pxx6;
% load power_Deltaband_d117Preeo 
% Pxx_PreNone9=Pxx6;
load power_Deltaband_d118Preeo 
Pxx_PreNone10=Pxx6;
% load power_Deltaband_d110Preeo 
% Pxx_PreNone11=Pxx6;
load power_Deltaband_d120Preeo 
Pxx_PreNone12=Pxx6;

Pxx_PreNone_Delta=[Pxx_PreNone1;Pxx_PreNone5;Pxx_PreNone8;Pxx_PreNone10;Pxx_PreNone12];
Pxx35=mean(Pxx_PreNone_Delta);

subplot(1,3,2)
topoplot(Pxx35,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Non responsive-Delta Band';'Before'})
h=colorbar;
h.Location = 'southoutside';
%%%%%%%%%%%%%%%%%%% Difference Power
Pxx208=Pxx35-Pxx36;
subplot(1,3,3)
topoplot(Pxx208,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Delta Band';'Difference Power'})
h=colorbar;
h.Location = 'southoutside';
%% Theta

load chanlocs19
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  Pre Reponsive

% load power_Thetaband_219Preeo 
% Pxx_PreNone1=Pxx7;
load power_Thetaband_312Preeo 
Pxx_PreNone2=Pxx7;
load power_Thetaband_313Preeo 
Pxx_PreNone3=Pxx7;
load power_Thetaband_112Preeo 
Pxx_PreNone4=Pxx7;
% load power_Thetaband_316Preeo 
% Pxx_PreNone5=Pxx7;
load power_Thetaband_326Preeo 
Pxx_PreNone6=Pxx7;
load power_Thetaband_115Preeo 
Pxx_PreNone7=Pxx7;
% load power_Thetaband_116Preeo 
% Pxx_PreNone8=Pxx7;
load power_Thetaband_117Preeo 
Pxx_PreNone9=Pxx7;
% load power_Thetaband_118Preeo 
% Pxx_PreNone10=Pxx7;
load power_Thetaband_110Preeo 
Pxx_PreNone11=Pxx7;
% load power_Thetaband_120Preeo 
% Pxx_PreNone12=Pxx7;

Pxx_PreNone_Theta=[Pxx_PreNone2;Pxx_PreNone3;Pxx_PreNone4;Pxx_PreNone6;Pxx_PreNone7;Pxx_PreNone9;Pxx_PreNone11];
Pxx38=mean(Pxx_PreNone_Theta);

subplot(1,3,1)
topoplot(Pxx38,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Responsive-Theta Band';'Before'})
h=colorbar;
h.Location = 'southoutside';
%%%%%%%%%%%% Non-responsive
load power_Thetaband_219Preeo 
Pxx_PreNone1=Pxx7;
% load power_Thetaband_312Preeo 
% Pxx_PreNone2=Pxx7;
% load power_Thetaband_313Preeo 
% Pxx_PreNone3=Pxx7;
% load power_Thetaband_112Preeo 
% Pxx_PreNone4=Pxx7;
load power_Thetaband_316Preeo 
Pxx_PreNone5=Pxx7;
% load power_Thetaband_326Preeo 
% Pxx_PreNone6=Pxx7;
% load power_Thetaband_115Preeo 
% Pxx_PreNone7=Pxx7;
load power_Thetaband_116Preeo 
Pxx_PreNone8=Pxx7;
% load power_Thetaband_117Preeo 
% Pxx_PreNone9=Pxx7;
load power_Thetaband_118Preeo 
Pxx_PreNone10=Pxx7;
% load power_Thetaband_110Preeo 
% Pxx_PreNone11=Pxx7;
load power_Thetaband_120Preeo 
Pxx_PreNone12=Pxx7;

Pxx_PreNone_Theta=[Pxx_PreNone1;Pxx_PreNone5;Pxx_PreNone8;Pxx_PreNone10;Pxx_PreNone12];
Pxx39=mean(Pxx_PreNone_Theta);

subplot(1,3,2)
topoplot(Pxx39,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Non Responsive-Theta Band';'Before'})
h=colorbar;
h.Location = 'southoutside';
%%%% Difference Power

Pxx210=Pxx39-Pxx38;
subplot(1,3,3)
topoplot(Pxx210,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Theta Band';'Difference Power'})
h=colorbar;
h.Location = 'southoutside';
%% Alpha
% Post
load chanlocs19
%Pre Responsive
% load power_Alphaband_219Preeo 
% Pxx_PreNone1=Pxx8;
load power_Alphaband_312Preeo 
Pxx_PreNone2=Pxx8;
load power_Alphaband_313Preeo 
Pxx_PreNone3=Pxx8;
load power_Alphaband_112Preeo 
Pxx_PreNone4=Pxx8;
% load power_Alphaband_316Preeo 
% Pxx_PreNone5=Pxx8;
load power_Alphaband_326Preeo 
Pxx_PreNone6=Pxx8;
load power_Alphaband_115Preeo 
Pxx_PreNone7=Pxx8;
% load power_Alphaband_116Preeo 
% Pxx_PreNone8=Pxx8;
load power_Alphaband_117Preeo 
Pxx_PreNone9=Pxx8;
% load power_Alphaband_118Preeo 
% Pxx_PreNone10=Pxx8;
load power_Alphaband_110Preeo 
Pxx_PreNone11=Pxx8;
% load power_Alphaband_120Preeo 
% Pxx_PreNone12=Pxx8;

Pxx_PreNone_Alpha=[Pxx_PreNone2;Pxx_PreNone3;Pxx_PreNone4;Pxx_PreNone6;Pxx_PreNone7;Pxx_PreNone9;Pxx_PreNone11];
Pxx40=mean(Pxx_PreNone_Alpha);

subplot(1,3,1)
topoplot(Pxx40,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Responsive-Alpha Band';'Before'})
h=colorbar;
h.Location = 'southoutside';
%%%%%%%% Non Responsive
load power_Alphaband_219Preeo 
Pxx_PreNone1=Pxx8;
% load power_Alphaband_312Preeo 
% Pxx_PreNone2=Pxx8;
% load power_Alphaband_313Preeo 
% Pxx_PreNone3=Pxx8;
% load power_Alphaband_112Preeo 
% Pxx_PreNone4=Pxx8;
load power_Alphaband_316Preeo 
Pxx_PreNone5=Pxx8;
% load power_Alphaband_326Preeo 
% Pxx_PreNone6=Pxx8;
% load power_Alphaband_115Preeo 
% Pxx_PreNone7=Pxx8;
load power_Alphaband_116Preeo 
Pxx_PreNone8=Pxx8;
% load power_Alphaband_117Preeo 
% Pxx_PreNone9=Pxx8;
load power_Alphaband_118Preeo 
Pxx_PreNone10=Pxx8;
% load power_Alphaband_110Preeo 
% Pxx_PreNone11=Pxx8;
load power_Alphaband_120Preeo 
Pxx_PreNone12=Pxx8;

Pxx_PreNone_Alpha=[Pxx_PreNone1;Pxx_PreNone5;Pxx_PreNone8;Pxx_PreNone10;Pxx_PreNone12];
Pxx41=mean(Pxx_PreNone_Alpha);

subplot(1,3,2)
topoplot(Pxx41,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Non Responsive-Alpha Band';'Before'})
h=colorbar;
h.Location = 'southoutside';

%Difference Power
Pxx220=Pxx41-Pxx40;
subplot(1,3,3)
topoplot(Pxx220,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Alpha Band';'Difference Power'})
h=colorbar;
h.Location = 'southoutside';
%% Beta
% Post
load chanlocs19
%Pre REsponsive
% load power_Betaband_219Preeo 
% Pxx_PreNone1=Pxx9;
load power_Betaband_312Preeo 
Pxx_PreNone2=Pxx9;
load power_Betaband_313Preeo 
Pxx_PreNone3=Pxx9;
load power_Betaband_112Preeo 
Pxx_PreNone4=Pxx9;
% load power_Betaband_316Preeo 
% Pxx_PreNone5=Pxx9;
load power_Betaband_326Preeo 
Pxx_PreNone6=Pxx9;
load power_Betaband_115Preeo 
Pxx_PreNone7=Pxx9;
% load power_Betaband_116Preeo 
% Pxx_PreNone8=Pxx9;
load power_Betaband_117Preeo 
Pxx_PreNone9=Pxx9;
% load power_Betaband_118Preeo 
% Pxx_PreNone10=Pxx9;
load power_Betaband_110Preeo 
Pxx_PreNone11=Pxx9;
% load power_Betaband_120Preeo 
% Pxx_PreNone12=Pxx9;

Pxx_PreNone_Beta=[Pxx_PreNone2;Pxx_PreNone3;Pxx_PreNone4;Pxx_PreNone6;Pxx_PreNone7;Pxx_PreNone9;Pxx_PreNone11];
Pxx42=mean(Pxx_PreNone_Beta);

subplot(1,3,1)
topoplot(Pxx42,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Responsive-Beta Band';'Before'})
h=colorbar;
h.Location = 'southoutside';
%%% non responsive
load power_Betaband_219Preeo 
Pxx_PreNone1=Pxx9;
% load power_Betaband_312Preeo 
% Pxx_PreNone2=Pxx9;
% load power_Betaband_313Preeo 
% Pxx_PreNone3=Pxx9;
% load power_Betaband_112Preeo 
% Pxx_PreNone4=Pxx9;
load power_Betaband_316Preeo 
Pxx_PreNone5=Pxx9;
% load power_Betaband_326Preeo 
% Pxx_PreNone6=Pxx9;
% load power_Betaband_115Preeo 
% Pxx_PreNone7=Pxx9;
load power_Betaband_116Preeo 
Pxx_PreNone8=Pxx9;
% load power_Betaband_117Preeo 
% Pxx_PreNone9=Pxx9;
load power_Betaband_118Preeo 
Pxx_PreNone10=Pxx9;
% load power_Betaband_110Preeo 
% Pxx_PreNone11=Pxx9;
load power_Betaband_120Preeo 
Pxx_PreNone12=Pxx9;

Pxx_PreNone_Beta=[Pxx_PreNone1;Pxx_PreNone5;Pxx_PreNone8;Pxx_PreNone10;Pxx_PreNone12];
Pxx41=mean(Pxx_PreNone_Beta);

subplot(1,3,2)
topoplot(Pxx41,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Non Responsive-Beta Band';'Before'})
h=colorbar;
h.Location = 'southoutside';

%%%Difference Power
Pxx222=Pxx41-Pxx42;
subplot(1,3,3)
topoplot(Pxx222,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Beta Band';'Difference Power'})
h=colorbar;
h.Location = 'southoutside';
%% Gamma
% Post
load chanlocs19
%Pre responsive
% load power_Gammaband_219Preeo 
% Pxx_PreNone1=Pxx10;
load power_Gammaband_312Preeo 
Pxx_PreNone2=Pxx10;
load power_Gammaband_313Preeo 
Pxx_PreNone3=Pxx10;
load power_Gammaband_112Preeo 
Pxx_PreNone4=Pxx10;
% load power_Gammaband_316Preeo 
% Pxx_PreNone5=Pxx10;
load power_Gammaband_326Preeo 
Pxx_PreNone6=Pxx10;
load power_Gammaband_115Preeo 
Pxx_PreNone7=Pxx10;
% load power_Gammaband_116Preeo 
% Pxx_PreNone8=Pxx10;
load power_Gammaband_117Preeo 
Pxx_PreNone9=Pxx10;
% load power_Gammaband_118Preeo 
% Pxx_PreNone10=Pxx10;
load power_Gammaband_110Preeo 
Pxx_PreNone11=Pxx10;
% load power_Gammaband_120Preeo 
% Pxx_PreNone12=Pxx10;

Pxx_PreNone_Gamma=[Pxx_PreNone2;Pxx_PreNone3;Pxx_PreNone4;Pxx_PreNone6;Pxx_PreNone7;Pxx_PreNone9;Pxx_PreNone11];
Pxx44=mean(Pxx_PreNone_Gamma);

subplot(1,3,1)
topoplot(Pxx44,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Responsive-Gamma Band';'Before'})
h=colorbar;
h.Location = 'southoutside'
%%%%% Pre Non reponsive
load power_Gammaband_219Preeo 
Pxx_PreNone1=Pxx10;
% load power_Gammaband_312Preeo 
% Pxx_PreNone2=Pxx10;
% load power_Gammaband_313Preeo 
% Pxx_PreNone3=Pxx10;
% load power_Gammaband_112Preeo 
% Pxx_PreNone4=Pxx10;
load power_Gammaband_316Preeo 
Pxx_PreNone5=Pxx10;
% load power_Gammaband_326Preeo 
% Pxx_PreNone6=Pxx10;
% load power_Gammaband_115Preeo 
% Pxx_PreNone7=Pxx10;
load power_Gammaband_116Preeo 
Pxx_PreNone8=Pxx10;
% load power_Gammaband_117Preeo 
% Pxx_PreNone9=Pxx10;
load power_Gammaband_118Preeo 
Pxx_PreNone10=Pxx10;
% load power_Gammaband_110Preeo 
% Pxx_PreNone11=Pxx10;
load power_Gammaband_120Preeo 
Pxx_PreNone12=Pxx10;

Pxx_PreNone_Gamma=[Pxx_PreNone1;Pxx_PreNone5;Pxx_PreNone8;Pxx_PreNone10;Pxx_PreNone12];
Pxx43=mean(Pxx_PreNone_Gamma);

subplot(1,3,2)
topoplot(Pxx43,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Non Responsive-Gamma Band';'Before'})
h=colorbar;
h.Location = 'southoutside';
%%%%%Difference Power
Pxx223=Pxx43-Pxx44;
subplot(1,3,3)
topoplot(Pxx223,EEG.chanlocs,'electrodes','labelpoint','style','both','colormap',colormap(jet(256)),'chaninfo',EEG.chaninfo)
title ({'\fontsize{16}Gamma Band';'Difference Power'})
h=colorbar;
h.Location = 'southoutside',